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The voltage- current (V-I) characteristics of tungsten (W) strips patterned by focused-ion-beam assisted 
deposition were studied at different temperatures and various magnetic fields applied perpendicular to the 
strips. At fields of 1 Tesla (T) and 2 T, a scaling behavior was observed in V-I isotherms, which can be fitted 
by quasi-three-dimensional (quasi-3D) vortex-glass (VG) to liquid transition theory. 

The existence of vortices in the mixed state of a superconductor is one of the most fascinating research topics 
in superconducting physics 1 " 16 . The vortices can form liquid, crystalline, and glassy phases. The discovery of 
high critical temperature (T c ) superconductors (HTS) has made the detailed study of vortex phases possible 
because of the large anisotropy and thermal fluctuations in HTS systems 1 " 7 . In the presence of strong disorder, the 
zero resistance state corresponds to a vortex glass (VG) phase and with increasing temperature the transition to 
the vortex liquid phase induces a dissipative state 2,17 . This transition has been demonstrated by transport mea- 
surement in HTS systems like YBa 2 Cu 3 0 7 (YBCO) 1 ' 3 , Pr 2 . x Ce x Cu0 4 . y (PCCO) 7 and Bi 2 Sr 2 . x La x Cu0 6 + s 
(BSLCO) 4 ' 6 . However, VG transition-like behavior in conventional low critical temperature superconductors 
(LTS) has only been observed in very few systems, like Nb 8 " 10 , In 11 , and Mo 3 Si 13 films. 

In the last few decades, a technique - focused ion beam induced deposition (IBID) has been developed and 
successfully used to fabricate low dimensional superconducting materials including W 18 " 23 . The reported T c of 
amorphous W-based nanostructures fabricated by IBID is around 5.2 K 18 , which is much higher than that of 
single-crystal bulk W (12 mK) 24 . This synthesis technique and the associated increase in T c makes it possible to 
observe different vortex states in mesoscopic W system. In two dimensional W systems fabricated by IBID, the 
direct observation of vortex melting utilizing scanning tunneling microscopy/spectroscopy (STM/STS) has been 
reported 25 . We report here a systematic study of V-I properties of W nanostructures. 

In this paper, W strips with different dimensions are fabricated using a focused ion beam- chemical vapor 
deposition (FIB-CVD) system. Experimental studies of electronic transport properties, especially the V-I curve 
measurement of individual W strips at various temperatures and magnetic fields are reported. Under zero 
magnetic field, the strip jumps from superconducting to normal state directly as the current is increased to a 
temperature dependent critical value. Under a finite magnetic field at temperatures below a specific value T g , the 
jump to the high resistance state is preceded by a state of low residual resistance. The size of the jump in the V-I 
curve decreases with increasing temperature and vanishes as T approaches T g . For T > T g , the V-I curves are 
continuous resembling a sigmoid function with a maximum curvature that decreases with increasing temper- 
ature. When plotted in log-log scale, the V-I curve at T = T g shows a power law dependence. Curves at T > T g show 
a concave curvature with increasing I and curves measured at T < T g prior to the jump show a convex curvature. 
This evolution in the shape of the V-I curves is similar to that found in HTS and type II superconducting films 
mentioned above, which has been explained by VG transition theory. 

Results 

A commercial FIB etching and deposition system (FIB/SEM FEI Quanta 200 3D) utilizing a beam of 30 keV 
singly- charged Ga + ions was used for sample synthesis. The substrate used was a lightly p-doped Si wafer with a 
1 um thick Si 3 N 4 insulating layer on the surface. Tungsten hexacarbonyl [W(CO) 6 ] gas was injected on the 
sample surface near the focused beam scanned area through a needle nozzle. The ion beam was then scanned over 
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the sample surface in the designed pattern, depositing the non-vol- 
atile component of the precursor on the surface. The volatile com- 
ponent was removed and the base pressure before introducing the 
W(CO) 6 gas was lower than 1 X 10~ 5 Torr. 

During deposition the system pressure was in the range of ~ 10" 5 
Torr. Lateral nanostrips with various lengths, widths and thicknesses 
were deposited by a 30 pA ion-beam current onto the substrate. 
Changing the ion-beam current or the volume per dose can change 
the superconducting properties of the resulting W strips. The sample 
W strips were connected to large Pt contact pads, which serve as 
electrodes for four-terminal transport measurements, with four W 
strips whose deposition parameters were the same as that of the 
sample. The normal Pt electrodes were also synthesized using the 
IBID technique. A scanning electron micrograph (SEM) of a W 
nanostrip with 4 connecting W strips can be seen in the inset of 
Figure 1(a). The volume per dose used in deposition is 6.8 X 
10" 1 unrVnC for samples 1 and 3, and 3.4 X 10" 1 um 3 /nC for sample 
2. The length, width and thickness of the samples were measured by a 
scanning electron microscope. Transport measurements were car- 
ried out in a physical property measurement system (PPMS- 
Quantum Design). The data from three samples are shown in this 
paper. The sample dimensions are 1.43 um long, 250 nm wide and 
50 nm thick (sample 1); 3.2 um long, 250 nm wide and 50 nm thick 
(sample 2) and 3.5 um long, 200 nm wide and 30 nm thick (sample 
3). The thickness of a similarly deposited W strip was measured by 
first cutting the strip with ion beam and then examined in a SEM 
chamber at a tilted angle. We estimate the uncertainty in the thick- 
ness to be no more than 5 nm. The elemental concentrations of our 
nanostrips detected by energy dispersive spectrum (EDS) are about 
40%, 40%, and 20%, for W, C, and Ga, respectively, which is con- 
sistent with the value detected by electron probe microanalysis 
(EPMA) 18 . 

An SEM image of sample 1 is shown in the inset of Figure 1(a). The 
length of the sample is defined as the distance L between the inner 
edges of the two voltage electrodes. Figure 1(a) shows resistance as a 
function of temperature (R-T) for sample 1 at an excitation current of 
50 nA in the temperature range of 2-300 K at zero magnetic field 
(H). This sample exhibits a superconducting transition and zero 
residual resistance below T c ~ 4.8 K. The sample does not show 
metallic behavior in the normal state and the resistance actually 
shows a 10% increase when it is cooled from 300 K to 8 K. This 
non-metallic behavior is likely a consequence of the presence of 
impurities such as carbon, Ga and oxides and also high density of 
defects in the amorphous W strips formed during the FIB-CVD 
fabrication process 18 ' 23 . In order to study the superconductivity near 
or below T c , the R-T curves in the low temperature regime (2 K ~ 
6 K) are measured. Figure 1(b) shows £-T curves in the low-temper- 
ature regime measured at different H. In this paper, the applied 
magnetic fields are always perpendicular to the plane of the sample 
and substrate. The R-T curve of the W strip at zero field exhibits a 
superconducting transition at T c = 4.9 K (onset T c ) and the width of 
the transition is about 0.2 K. On increasing H, the superconducting 
transition broadens and shifts to lower temperatures. Figure 1(c) 
shows the resistance of the W strip as a function of H at different 
temperatures. A superconducting transition with an onset critical 
field (H C2 ) larger than 8 T at 2 K is seen. The H C2 of the strip 
decreases with increasing temperature. At a temperature of 6.0 K, 
the superconducting transition is barely noticeable. 

Discussion 

Although bulk W is a type I superconductor, the W strips measured 
in these experiments can be thought of as a type II conventional 
superconductor due to its finite thickness and amorphous structure 
with high density of defects as seen in transmission electron micro- 
scope (TEM) studies 23 . The disorders and defects act as pinning 
centers for vortices. For type II conventional superconductor 
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Figure 1 | Transport measurements of the W strip, (a) Resistance vs. 
temperature of W strip (sample 1). Inset is the SEM image of the sample 
with 4 connecting W electrodes; (b) Resistance vs. temperature of the W 
strip in different magnetic fields applied perpendicular to the plane of the 
sample; (c) Magnetoresistance of the sample at different temperatures. 



Nb 8 " 10 , the electric field - current density (E-J) isotherms measured 
at various temperatures at a given field value show qualitatively dif- 
ferent behavior demarcated by a VG transition temperature, T g . Near 
T ff the curves show a power law relationship. For T > Tg, the E-J 



curves exhibit ohmic behavior and the resistivity remains constant at 
low currents. For T < Tg, E ~ exp (-J/J^ 1 ) resulting in a convex curv- 
ature in the E-J curves which is a signature of the vortex glass phase 6 . 

According to the VG transition theory 6 , the V-I curves at different 
temperatures near T g can be scaled into two different branches by the 
scaling law, 
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y//(r-r,)" (z+2 - D) =A ± (//|r-r,| 0(D - 1) ) (i) 

where v is the exponent of the coherence length £~ \ T— T g \ °, 
and z is the dynamical exponent of the coherent time ^respectively, 
D is the dimension of the system, and A + are the scaling functions 
above and below T g . In the vortex liquid phase for the temperature 
above T g , the linear resistivity p z/n at low current behaves as, 

ft „*(r-r,r (z+2 - D) (2) 

The V-I isotherm at T g should satisfy the power law, 

V(I)| r _^xI a (3) 

with a=(z+l)/(D— 1). Previous theoretical and experimental 
studies show that z = 4 ~ 7 and u = 1 ~ 2 are reasonable values for a 



3D VG phase transition 1,6,7,9 . In recent years, quasi-2D VG transi- 
tions have been reported, and the ranges of z and v are similar to 
those for 3D 12 . 

The V-I curves of the W strips are obtained at various magnetic 
fields and different temperatures. In this paper, we have measured 
the V-I curves by ramping the current from low to high and from 
high to low values. Identical results are found and we report only the 
curves collected when the current is being scanned downward. 
Figures 2 (a) -2(c) show V-I curves measured between 2 and 6 K at 
applied fields of 0 Oe, 10 kOe and 20 kOe respectively. At zero field, 
the measured voltage jumps abruptly from zero to the normal state 
value at well-defined temperature dependent critical current. 
However, at a field of 10 kOe (Figure 2(b)), at temperature below 
Tg, V first increases slowly with I. This is followed by an abrupt jump 
and then a more gradual increase with further increase in I towards 
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Figure 2 | V-I measurement of the W strip at different temperature and under magnetic field, (a-c) V-I characteristics measured at temperatures 
ranging from 2 to 6 KforH = 0 kOe y 10 kOeand20 kOe y respectively. The red arrow points out that the measurement temperature changes from 2.0 Kto 
6.0 K. (d-f) V-I characteristics in log-log scales The dashed lines in panels (e) and (f) highlight power law V-I dependence at T = T g (4.05 K for 
measurements taken at 10 kOe and 3.4 K at 20 kOe). 
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ohmic behavior. V-I curves measured at T > T g do not show an 
abrupt jump and the curves are continuous resembling a sigmoid 
function. The maximum curvature of these curves decreases with 
increasing temperature. Similar evolution of V as a function of the 
current is seen at a field of 20 kOe (Figure 2(c)). The value of T g 
decreases with magnetic field. The V-I curves shown in Figures 2(b) 
and 2(c) show that T g at 10 kOe is near 4.0 K and at 20 kOe is near 
3.4 K. We will show below that T g can be more precisely identified by 
a more quantitative procedure. We define the critical current I C o> 
where the voltage changes from zero to non-zero and critical current 
I Ci , where the voltage jumps rapidly. 

The V vs. I dependence at low current is revealed in more details in 
log-log plots. Figure 2(d) clearly shows that at zero field the V-I curves 
show an abrupt change from the superconducting to the normal 
ohmic behavior at well defined I c s for all T < T c . For measurements 
at fields of 10 and 20 kOe, the log V - log I curves (Figures 2(e) and 
2(f)) at currents smaller than I C1 show a smooth evolution in the 
curvature from being convex to concave with increasing temper- 
ature. The curves with different curvature are separated by a curve 
showing power law dependence measured at or near T = T g . This 
behavior is reminiscent of the vortex-glass state with the crossover 
temperature T C _ Q corresponding to the VG transition temperature - 
T g . When the temperature is below T g , and the sample is in the vortex 
glass state, the resistance falls rapidly with decreasing current and is 
zero below I co . At temperatures above T g , the sample is in the vortex 
liquid state, the V-I curves exhibit ohmic behavior and the resistance 
remains constant even at small excitation currents. T a can be more 




Figure 3 | dV/dl-I plots at temperatures ranging from 2 to 6 K for 10 kOe 
and 20 kOe, respectively. The insets show the inverse of the maximum 
values of dV/dl (black line) and the changing rate of the slope of the black 
line. The blue line exhibits that the max changing rate is at 4.05 K and 3.4 K 
for 10 kOe and 20 kOe, respectively. 



precisely identified by plotting the slope (dV/dl) as a function of I for 
each T numerically computed from the data of Figures 2(b) and 2(c). 
For T < T g , the curves show a sharp and theoretically infinite dV/dl 
peak at some current, for T > T g , the dV/dl shows a low and wide 
peak (shown in Figures 3(a) and 3(b)). After plotting the inverse of 
the maximum dV/dl values shown in insets of Figures 3(a) and 3(b) 
for each temperature (black line), the changing rate of the slope can 
be calculated as a function of temperature shown in blue line in the 
inset of Figures 3(a) and 3(b). The T g of sample is the point which has 
the maximum changing rate of inversed maximum dV/dl value vs. T 
(4.05 K for the curves under a field of 10 kOe and 3.4 K for 20 kOe). 
A red dashed line in Figures 2(e) and 2(f) represents the V-I behavior 
at T g . The fitted line satisfies the relation — V ~ I a where a is 3.15 and 
3.16 for Figures 2(e and 2(f)) respectively. The values at different 
fields are in close agreement. 

The inset of Figure 4 shows the logR-T curve of sample 1 at 10 kOe 
with fit based on Eq. (2). The T g value is 4.05 K. The slope of the 
fitting curve is o(z + 2 — D) . In order to get reasonable values of z and 
v within the limits mentioned above (z = 4~7 and u=l~2), the 
dimension D is taken to be 3 and the fitting parameter values are 
found to be z = 5.31 and v = 0.99. This means our W strip system 
shows quasi- 3D behavior. This is reasonable since the coherence 
length of W strip is about 6 nm 17 , which is much smaller than the 
thickness and width of the sample. Figure 4(a) shows the scaling 
collapse of the V-I isothermals according to Eq. (1) based on the 
experimental V-I curves at 10 kOe and the obtained parameter values 
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Figure 4 | Quasi- 3 D VG scaling of the V-/ curves measured at 10 icOeand 
20 kOe respectively. The insets show a logarithmic plot of the 
temperature dependence of the linear resistivity along with the fit used to 
obtain the scaling parameters. 
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(z = 5.3 land v = 0.99). The experimental V-I curves formed two uni- 
versal branches corresponding to the data above and below T g 
(4.05 K). Similar universal scaling behavior can be obtained for the 
V-I curves at a field of 20 kOe (Figure 4(b)) with parameter values of 
T g = 3AK, z = 5.33 and v = 0.99. 

The V-I isothermals under different magnetic fields at 2 K 
(Figure 5(a)) are qualitatively similar to the V-I curves at different 
temperatures at fixed magnetic fields. The difference is that the V-I 
curves are separated into two parts by a particular field - H g shown by 
a red dashed line in Figure 5(a). The isotherms are concave at the 
fields higher than H g , while they are convex at fields lower than H g . 
H g plays the same role for these curves as T g did in the discussion 
above. When the applied field is smaller than H g , the voltage drops 
rapidly with decreasing current. At fields larger than H g , the V-I 
curves exhibit ohmic behavior and the resistance remains constant 
in small current. This behavior exhibits the characteristics of a vortex 
glass-liquid transition. H g can be considered to be the critical field of 
the VG phase transition. On the basis of the VG transition theory, the 
Eq. (1) can be rewritten in the following form: 



V I H-H 0 



xuO + 2-D) 



=a ± (//|h-h/ d " 1) ) 



(4) 



Utilizing Eq. (4) and an approach similar to the one described above 
(also see Supplementary Information), the related parameter values 
can be obtained (H g = 35 kOe, z = 7.07 and v = 0.89), the V-I iso- 
thermals collapse onto two different branches as shown in 
Figure 5(b). The values of z and v are a little deviation from the 
reasonable ranges mentioned above, which is possibly due to the 
limited thickness and width of the nanostrip. 



Figure 6(a) is the H-T phase diagram for sample 1 (W strip). The 
values of the upper critical field H C2 , are extracted from the R-H 
curves at different temperatures shown in Figure 1(c). The fitting 
line in Figure 6(a) is based on the empirical temperature dependence 
ofH c (Eq.(5)). 

Jc) 



H C 2(T)=H C2 (0) 



1- 



(5) 



In Figure 6(a), the black dots are T g at 1 T and 2 T obtained from the 
V-I curves in Figure 3. The dashed fitting line is obtained using 
empirical formula (6) 17 , 

H m (T)~(l-^Q 2 " (6) 

H m (T) is the field corresponding to the vortex glass melting temper- 
ature (T g ). From the V-I data at different temperatures and under 
different magnetic fields (see Figures 2(e) & 2(d), Figure 5(a)), the 
value of T g at 10 kOe and 20 kOe, H g dX2 K can be derived. With two 
values of parameter v obtained by fitting Eq. (2) earlier (0.99 and 0.89 
under 10 kOe and 20 kOe field respectively), the fitting lines with two 
v values are almost the same. We found that the data point -H g dX2 K 
(the star mark) in Figure 6(a) lies on the red dashed fitting line which 
is the VG melting line. According to the phase diagram, each T g 
corresponds a unique field along the transition curve. Thus, the H g 
should be equivalent to T g on the melting line. 

According to Figure 2, the W strip is in totally superconducting 
state (zero voltage) below J co , then the voltage increases slowly till I C1 
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Figure 5 | V-I curves and Quasi-3D VG of the W strip at 2 K. (a) V-I 

characteristics measured at 2 K for magnetic fields ranging from 5 to 
80 kOe; (b) Quasi-3D VG scaling of the V-I curves measured at 2 K. 
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phase diagram of magnetic field dependence of the melting temperature, 
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symbols. For 0 Oe data, the values of J co and J C1 are the same. 
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Figure 7 | Quasi-3D VG scaling of the V-I curves measured at 2 K for 
sample 2 and 3. (a) Quasi-3D VG scaling of the V-I curves measured at 2 K 
for sample 2; (b) Quasi-3D VG scaling of the V-I curves measured at 2 K 
for sample 3. 

and jumps to the normal state. Figure 6(b) shows the critical current 
density J co and J Ci , calculated from I co and I CI respectively, changing 
with temperature at 0, 10 and 20 kOe. The J co vs T curves obey the 
empirically parabolic law (Eq. (7)). For 0 Oe, the J co and J C1 are the 
same. The J C o is 3.2 x 10 5 A / cm 2 at 2 K under 0 Oe field (The value in 
P - W film, of which T c is 3.35 K, is about 10 6 - 10 7 A/ cm 2 26 ). With 
the temperature and magnetic field increasing, the values of J c 
decreases. 



/ Ci (r)=/ Ci (o) 



l- 



(7) 



The results discussed above were verified on the other two W strip 
samples by measuring their V-I characteristics at 2 K at different 
applied H. The crossover behavior (the curvature changes from con- 
vex to concave) discussed above was also observed in these samples at 
a particular field (H g ). In Figures 7(a) and 7(b), the V-I isothermals 
for samples 2 and 3 are scaling collapsed onto two branches obtained 
by using parameters: H g = 30.5 kOe, z = 5.24 and u = 0.67 for sample 
2, H g = 28 kOe, z = 6.02 and v = 0.53 for sample 3. According to 
Figure 7, sample 2 obeys the VG scaling better than sample 3. This 
may be because that the thickness and width of sample 3 are dis- 
tinctly smaller than those of sample 1 and sample 2. The small thick- 
ness and width limit the applicability of quasi- 3D and quasi-2D VG 
transition theory to sample 3 and quasi- ID size effects may need to 
be taken into account in this situation. 

In summary, the V-I curves of the W strips fabricated by focused 
ion beam & chemical vapor deposition are measured at different 



magnetic fields and temperatures. The scaling analysis of the V-I 
isothermals of W nanostrip (sample 1 and sample 2) follows the 
quasi- 3D VG transition theory well. Every VG transition temper- 
ature T g has a corresponding H g along the VG melting line i.e. H g 
appears to be equivalent to T g in VG transition. However, the thin- 
nest W strip (sample 3) shows a larger deviation from the conven- 
tional quasi-3D VG fitting. This may be because quasi- ID effects are 
not negligible in this situation 27 and have to be considered in future 
theoretical studies. 
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